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We report on a study of the finite-temperature chiral transition on an Nt — 4 lattice for 2 + 1 flavors of Kogut- 
Susskind quarks. We find the point of physical quark masses to lie in the region of crossover, in agreement with 
results of previous studies. Results of a detailed examination of the m u ,d = rn a case indicate vanishing of the 
screening mass of a meson at the end point of the first-order transition. 



1. Introduction 

An important issue in finite-temperature lat- 
tice QCD is the determination of the nature of 
the chiral phase transition for a realistic spec- 
trum of light up and down quarks and a heav- 
ier strange quark. Despite its importance, past 
studies of this "2+1" case have been few. For 
the Kogut-Susskind action, all of them have been 
made around 1990 1§-§|. 

It was found in these studies that the chiral 
phase transition changes from a first-order tran- 
sition to a crossover as the strange quark mass m s 
increases beyond a critical value m c s for a fixed de- 
generate up and down quark mass m Ut d, in agree- 
ment with predictions of an effective a model of 
QCDj§[|. Results were also obtained |||§ which 
indicate the physical point of quark masses to 
lie in the region of crossover on the {m u ,d,rn s ) 
plane. However, these results were based on sim- 
ulations made at only a few sets of quark masses. 
Clearly a more extensive study is called for to 
have a full understanding of the phase structure 
in the 2-parameter space of (m u ^, m s ). Here we 
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report first results from our recent effort toward 
this goal. 

An interesting suggestion from a a model anal- 
ysis]^] is that the second-order transition ex- 
pected at the critical strange quark mass m c s is 
in the Ising universality class, with the massless 
mode provided by the a meson. A novel feature 
of our work is a study of the screening mass of a 
to examine this point. 

Our simulations are performed for the tempo- 
ral lattice size Nt = 4. An 8 3 x 4 lattice is em- 
ployed to make a survey of the phase structure 
varying j3, m Ut d and m s . A detailed investigation 
is then made along the flavor SU(3) symmetric 
line m Ut d — m s by another series of simulations 
on a 16 x 8 2 x 4 and a 16 3 x 4 lattice. For each pa- 
rameter set, (1-2) xlO 3 trajectories of unit length 
are generated by the hybrid R algorithm. 

2. Phase diagram on the (m U) d,m s ) plane 

We show the result of our phase diagram anal- 
ysis on an 8 3 x 4 lattice in Fig. [|. At m Ut d = 
m s = 0.01 a clear two-state signal is obtained by 
a comparison of runs with a hot and a cold start. 
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Figure 1. Phase diagram on (rn Ut d,m s ) plane ob- 
tained on a 8 3 x 4 lattice. Open symbols are 
results of Ref. ||] on a 16 3 x 4 lattice. Shaded 
line represents an estimate of the critical line, and 
cross the physical point (see text). 

At triangle points, we find a very sharp change of 
observables over a narrow range of /3, suggestive 
of a first-order transition, while only a smooth 
crossover is seen at diamond points. 

In Fig. ^ we plot results for the chiral con- 
densate along the line m Ut d = m s = to on a 
16 3 x 4 lattice (m < 0.03) or a 16 x 8 2 x 4 lat- 
tice (to > 0.04). On these lattices a two-state 
signal is found, which is clear at m = 0.01 and 
0.025, but less so at to — 0.03. The behavior 
above m = 0.04, on the other hand, indicates a 
crossover. 

The results taken together suggest that the 
critical line marking the end point of the first- 
order transition runs above the triangle points at 
mu,d < 0.025, and below the point m Uy d = m s = 
0.04 in Fig. |l|, which is indicated by the shaded 
line. If we assume T c — 150 MeV for the critical 
temperature, the physical point is located at the 
cross in Fig. [|, which is in the region of crossover. 
Our results are in agreement with those of the 
Columbia group obtained at m u> d = 0.025 on a 
16 3 x 4 lattice || . 

3. Results along the m u ,d = m s line 

We now discuss results along the line m u ,d — 
m s = m. Given our observation of a two-state 
signal at to — 0.01,0.025, one way to estimate 
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Figure 2. (ipip) as a function of (3 for various 
quark masses along the line m Ul d = m s = m. 



the value of the critical mass m c is to extrapolate 
the gap of the chiral condensate A (ipip) toward 
larger m where it vanishes. Employing the form 
A(ipip) oc (to c — to) 1 / 2 predicted by the mean- 
field analysis of the a model, we find m c ~ 0.034. 
If we use a naive linear extrapolation, we obtain 
to c ~ 0.049. A similar value m c w 0.045 was 
previously reported M by a linear extrapolation 
applied to old results|g||. 

In the region of crossover to > m c , we expect 
the peak height of the chiral susceptibility Xm to 
develop a singular behavior Xm oc (to — m c )~ z as 
to — ► to c . We calculate x™. f° r m > 0.04 with the 
histogram reweighting method. Assuming m c = 
0.034, we fit the peak height to the form above. A 
reasonable fit with x 2 /g?/ = 1.05 is obtained with 
the value of the exponent z = 0.67(3), which is 
comparable to the Ising value z ~ 0.79 and the 
mean- field value 2/3. 

In order to examine the screening mass M a of 
a meson, we employ 17(1) random source and no 
gauge fixing to evaluate the two quark loop contri- 
bution of the a propagator. Good results are ob- 
tained for the full a propagator with this method 
as illustrated in Fig. [|. 

The quark mass dependence of 7r and a screen- 
ing masses for m < 0.03, where we find a first- 
order transition, is plotted in Fig. f|(a). We ob- 
serve that M 2 decreases toward zero as to in- 
creases toward the critical value, both on the con- 
fined and the deconfined side of the transition, in 
contrast to M z which increases. 
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Figure 3. tt and a propagators for to — 0.1 at the 
transition point /3=5.25. 



Assuming a linear quark mass dependence, 
M% oc m c — to, predicted by the mean-field anal- 
ysis of the a model, we obtain m c — 0.034(3) in 
the confining phase and m c = 0.031(3) in the de- 
confining phase. These values are consistent with 
each other, and are also in agreement with the 
estimate from a square root extrapolation of the 
gap of (t^ip) discussed above. These results indi- 
cate vanishing of the a mass at the critical quark 
mass as suggested by the a model 

Results for larger quark masses (to > 0.04), 
where a crossover behavior is seen, is shown in 
Fig. ||(b). While M% decreases toward smaller 
values of to, the variation is too mild to attempt 
an independent estimate of m c . An interesting 
point which requires clarification is that M% stays 
considerably small compared to M% even for large 
quark masses. 



4. Conclusions and future work 

Our study with the Kogut-Susskind action sup- 
ports the previous conclusion with this action 
that there is no finite-temperature phase tran- 
sition for three flavors of quarks with physical 
masses. This means that a discrepancy with the 
conclusion from the Wilson action still remains. 

We also find a strong indication that the screen- 
ing mass of er vanishes at the end point of the 
first-order transition along the line m u .d = m s . 

We plan to extend analyses carried out here to 
the m Ut d 7^ m s case to further explore the real- 
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Figure 4. Screening masses M% and M% as a func- 
tion of to for (a) m < 0.03 and (b) to > 0.04 at 
the transition point. In (a) results on the two 
sides of first-order transition are shown. 

world QCD chiral transition. 
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